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GENERATORS AND ACCELERATORS 
By Gustav A. Carlson 
Ames Research Center 


SUMMARY 


Boundary layers on the walls of linear crossed-field plasma generators 
and accelerators cause circulating electrical currents that can greatly affect 
the operation of these devices. Expressions are derived for the velocity and 
current distributions in the open circuit device from the analysis for fully 
developed incompressible magnetohydrodynamic (MHD) flow in a square duct with 
nonconducting walls. From these distributions the viscous and ohmic heat gen- 
eration rates in the working fluid were calculated for Hartmann numbers up to 
24, where the magnetic force is dominant. Measured heat-transfer rates to 
the walls of an open circuit J * B plasma accelerator were found to increase 
with applied magnetic field strength at the insulator walls but to remain 
nearly constant at the electrode walls. A comparison of theory and experiment 
indicates that the calculated heat generation has both the magnitude and spa- 
tial distribution to explain the observed heat-transfer rates. The circula- 
ting currents significantly decelerate the flow near the center of the duct, 
and will exist in any MHD duct in which there is an appreciable boundary layer 
and a large conductivity extending to near the walls. 


INTRODUCTION 


In most cases the simple one-dimensional theory of magnetohydrodynamic 
duct flow overpredicts the actual performance of linear crossed-field plasma 
generators and accelerators. Most attempts to extend the theory so it will 
better predict actual performance have been directed toward the description of 
the current distribution in the plane perpendicular to the applied magnetic 
field. Factors which influence the current distribution in that plane include 
electrode segmentation, velocity and thermal boundary layers on cold elec- 
trodes, Hall effects, unequal electron and gas temperatures, and nonequilib- 
rium ionization. These factors have been shown to prevent the desired axial 
orientation of the J x B force (refs. 1-5). For open circuit operation of a 
J x B device, these investigations predict that the only magnetic effect will 
be an induced voltage proportional to the product of the applied magnetic 
field and the fluid velocity. Experimental measurements have departed from 
this prediction: The induced open circuit voltage is less than calculated 

from theory (refs. 6, 7); the impact pressure decreases with increasing mag- 
netic field; and the heat transfer to the insulator walls is greater than that 
estimated for ordinary boundary layers (ref. 8) . 


This degradation of performance is caused by electrical currents which 
circulate within the working fluid. Circulating current loops in a linear 
J x B device may extend in the axial direction (end looping), or may lie in 
the plane perpendicular to the direction of flow. End looping (refs. 9, 10) 
is caused by axial variations in the applied magnetic field and will not be 
treated in this paper. Circulating currents in the plane perpendicular to the 
direction of flow are caused by the nonuniform axial velocity distribution. 
Schneider and Wilhelm (ref. 7) utilize the well known Hartmann problem to show 
qualitatively how these circulating currents can cause induced open circuit 
voltage to have a nonlinear dependence on magnetic field. Reilly and Oates 
(ref. 8) discuss these circulating currents as a source of ohmic heating and 
as the cause of increased shear at the insulator walls; however, they report 
experimental measurements of local heat transfer for only one value of mag- 
netic field and present no calculations for boundary- layer heat transfer that 
include the effects of the circulating currents. 

The purposes of this paper are (1) to report the results of calculations 
of the heat generation in the working fluid of an open circuit J x B device 
that include the effects of the circulating currents, and (2) to compare these 
calculations to experimental measurements of bulk heat transfer to the duct 
walls for various values of applied magnetic field. It must be noted that the 
calculations are for heat generation, whereas the measurements are for heat 
transfer, and that an increase in heat generation exceeds the corresponding 
increase in heat transfer by the amount of thermal energy increase of the 
fluid. 


SYMBOLS 


a half-width of square duct 

B magnetic induction vector 

Bo applied magnetic induction 

E electric field strength vector 

H magnetic field strength vector 

H Hartmann number, B 0 a 

H x axial component of magnetic field strength 

3 current density vector 

jy current density in y direction 

j current density in z direction 

l duct length 
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static pressure 

rate of ohmic heat generation per unit volume 

rate of viscous heat generation per unit volume 

rate of heat generation or heat transfer 

velocity vector 

axial velocity 

average axial velocity 

Cartesian coordinates 

dimensionless coordinate, — 

a 

dimensionless coordinate, 

a 

dynamic viscosity of fluid 
density of fluid 

electrical conductivity of fluid 
dimensionless quantity 

THEORETICAL ANALYSIS 


The duct configuration and coordinate system for the circulating current 
model are shown in figure 1. In this model the applied magnetic field B 0 is 
assumed constant with only a z component. The axial velocity of the electri- 
cally conducting working fluid varies in both the y and z directions because 
of the boundary layers on the walls of the duct. This variation results in a 
nonzero value for curl (0 x B) ; however. Maxwell's electric field equation 
requires curl (E) = 0. It follows from Ohm's law,' J = a(E + 0 x B) , that 5 is 
not zero, even for the open circuit case. 

The circulating currents in the y-z plane of the open circuit duct are 
calculated from the analysis for fully developed incompressible MHD flow in a 
rectangular duct with nonconducting walls. The circulating current loops are 
assumed to remain entirely within the fluid as shown in figure 1. Such a cur- 
rent loop is more likely than one that includes current flow through the elec- 
trodes in the z direction (as assumed in the Hartmann problem) because the 
electrode-plasma interface resistance is appreciable. The nonconducting wall 
problem was solved exactly by Shercliff (ref. 11) in 1953 and has since been 
considered by others. For this paper, the exact solution presented by Hughes 
and Young (ref. 12) is the most useful. 
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Insulator wall 


Hughes and Young write the Navier-Stokes equation for fully developed, 
incompressible steady flow and combine the Maxwell equations V x E = 0 and 
V x H = J with Ohm's law for constant conductivity J = a(E + 0 x §) to form 
the steady-state magnetic transport equation. For the coordinate system of 
figure 1 the two final equations are 
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For flow in a rectangular duct with nonconducting walls, Hughes and Young show 
that the solution to equations (1) and (2) may be expressed as simple Fourier 
series. The dimensionless form of their solution is not well suited for calcu- 
lations of magnetic effect with flow rate held constant. With a choice of 
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dimensionless parameters better suited for such calculations, their solution 
for velocity, current, and pressure gradient in a square duct becomes 
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The local heat generation rate in the flow is divided into two parts, the 
viscous term q y = y[(3u/8y) 2 + (3u/3z) 2 ] and the ohmic term 

q. = (l/a)(jy 2 + j z 2 ) • These terms are calculated in dimensionless form from 
equations (3) through (5) assuming both a and y constant: 
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The total heat generation rates are found by integrating equations (7) and (8) 
over the volume of the duct. Applying Green's theorem and substituting from 
equations (1) and (2) simplifies the integration to yield the following solu- 
tion for total viscous and ohmic heat generation: 
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The integral in equations (9) and (10) is evaluated numerically by a double 
application of Simpson T s rule. The total heat generation rate (viscous plus 
ohmic) may be calculated without numerical integration, since the sum of 
equations (9) and (10) is simply 


* 

Q TO t 




(ii) 


The calculations for local and total heat generation rates were carried out 
for 11 values of Hartmann number, from 0.5 where the magnetic body force is 
negligible to 24 where it is dominant. 


EXPERIMENT 


The experiment was conducted on a linear J x B accelerator depicted 
schematically in figure 1. The accelerator channel is square at the entrance 
and diverges at a 2° half-angle in the x-y plane. There are 22 pairs of cop- 
per electrodes separated in the flow direction by 0.07-cm air gaps. The insu- 
lator walls are constructed of anodized aluminum and are insulated from the 
electrodes. As a further safeguard against shorting, each insulator wall of 
the channel consists of three slabs, with segmentation in the y direction. 

The electrodes and insulator walls are water cooled. The accelerator exhausts 
into a vacuum tank (operating pressure 2 . 5 mm Hg) and is fed by a supersonic 

constricted-arc jet (ref. 13] 1.27 cm in diameter and 30 cm long. The working 

fluid is unseeded argon added at the upstream end of the arc jet at a rate of 

3.65 g/s. The steady-state arc jet operating conditions are 600 A dc at 400 V. 

In the work reported here the accelerator was operated in the open circuit 
condition with values of applied magnetic field from 0 to 0.5 tesla (T) . 


6 


During each run a pitot probe and a stagnation heat-transfer rate probe 
were traversed in the y direction through the center of the plasma stream 
just downstream of the accelerator exit. Data from these probes were used to 
calculate the centerline velocity of the flow (refs. 14, 15). The velocity 
was independently checked by time of flight measurements of luminosity 
perturbations at the centerline of the plasma. 

The heat-transfer rate to each of the accelerator walls was recorded 
continuously by a galvanometer-type oscillograph. This was accomplished by 
delta-T transducers (ref. 16) which measured the temperature rise of the water 
in each coolant circuit (see fig. 2). The pressure drop through each coolant 
circuit was maintained constant from run to run, thus insuring an invariant 
flow rate. The calibration constants for the transducer-oscillograph system 
were obtained directly from known amounts of electrical power dissipated in 
sections of high electrical resistance tubing that form a permanent part of 
each coolant circuit. 


Electrode or 
insulator wall 



Signal to 
oscillograph 


Figure 2.- Typical coolant circuit. 

The measured wall heat -transfer rates were made dimensionless functions 
of the Hartmann number by use of the duct dimensions, a dynamic viscosity of 
2*10~ 4 N-s/m 2 , an electrical conductivity of 700 mho/m, and an average velocity 
of 2750 m/s. 


Q* = — (12) 
4u 2 y 2 


H = B 0 av^A7 (13) 

The transport properties were taken from De Voto (ref. 17) for equilibrium 
argon at 10 mm Hg pressure and 8000° K temperature, the approximate operating 
conditions in the accelerator. From the velocity relation for fully developed 
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flow in a square duct, u(y = z = 0) ^ 2 . 1 u, the average velocity was 
estimated as half the measured centerline velocity at zero magnetic field con- 
ditions. The fully developed flow relation was used because the traversing 
pitot probe measured a parabolic impact pressure profile, indicating that the 
wall boundary layers filled the duct. Although the quantitative comparison of 
experiment and theory is only as good as the estimates for viscosity, conduc- 
tivity, and average velocity, the qualitative comments which follow are not 
dependent on these estimates. 


COMPARISON OF THEORY AND EXPERIMENT 


The spatial distributions of the calculated ohmic and viscous heat 
generation rates for Hartmann numbers of 1 and 10 are shown in figures 3 and 4. 
At the smaller Hartmann number, the viscous heat generation distribution is 
nearly symmetric about the duct diagonal, and the ohmic generation is very 
small. At the larger Hartmann number, both the viscous and the ohmic heat gen- 
eration are much increased. (Note the change of scale between figs. 3 and 4.) 
These increases are caused by large circulating currents which develop as the 
Hartmann number is increased, with the largest current density appearing adja- 
cent to the insulator walls. These currents result in correspondingly large 
ohmic heat generation and, in addition, cause a flattening of the velocity dis- 
tribution in the z direction, which increases viscous heat generation 
because of the increased velocity gradient at the wall. 

The measured values of heat transfer at H = 0 were 6.79 kW to the insula- 
tor walls and 2.33 kW to each of the electrode walls. As H was increased to 
23, the measured heat transfer to the insulator walls increased to 25.2 kW, 
but the heat to the electrodes remained nearly constant. This concentration 
on the insulator walls of the heat-transfer increase with H is in excellent 
qualitative agreement with the calculated heat generation distribution 
depicted in figure 4. 

The calculated values for the integrated ohmic and viscous heat 
generation in the fluid and the experimentally measured heat transfer to the 
walls of the accelerator are compared in figure 5 as functions of the Hartmann 
number. To give all the functions a zero value at H = 0, the calculated vis- 
cous heat generation for H = 0 was subtracted from all calculated values of 
viscous heating, and the measured heat-transfer values for H = 0 were sub- 
tracted from all measured values of heat transfer. As expected, the measured 
total heat transfer increases less rapidly with H than does the calculated 
total heat generation (solid curve), because some of the heat generated is not 
transferred to the walls but remains in the fluid. The apparent leveling off 
of the measured heat transfer at the larger Hartmann numbers is not understood; 
additional experimentation with stronger magnetic fields will be required to 
determine whether the apparent trend persists. 

For nonzero values of the Hartmann number, the measured heat transfer to 
the upper electrode wall was consistently less than that to the lower elec- 
trode wall and, in some cases, was actually less than the heat transfer at 
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Figure 5.- Measured heat transfer and theoretical heat generation as a function of 

Hartmann number. 


H = 0. This behavior is probably due to the Hall effect, which causes the cir- 
culating current density near the center of the duct (where jy is negative) 
to have a component in the positive x direction. The resulting J x B force 
acts on the fluid in the negative y direction, tending to increase the pres- 
sure and the heat-transfer rate at the lower electrode wall. 

The circulating currents work to decelerate the supersonic flow near the 
center of the duct by two mechanisms: enthalpy increase and J x B decelera- 
tion. For H = 10, the calculated value for jy* at y = z = 0 is -1.42, 
which, for the approximate operating conditions of the experiment corresponds 
to a current density of 5.74 A/ cm 2 . The one-dimensional theory of MHD duct 
flow predicts that a current density of this magnitude and a magnetic field of 
0.214 tesla (H = 10) would result in a velocity decrease of 86 m/s for each 
centimeter of duct length. The measured and theoretical values of dynamic 
pressure at the stream centerline are shown in figure 6 as functions of the 
Hartmann number. The measured values were computed from the pitot probe data 
by means of the supersonic pitot tube formula, and the theoretical curve was 
calculated from equation (3) for y = z = 0. The theoretical decrease of 
dynamic pressure with increasing Hartmann number is in excellent agreement 
with the experimental measurements. 
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Figure 6.- Centerline dynamic pressure at the accelerator exit as a function of 

Hartmann number. 

Both for the theory and for reducing the measured heat-transfer data to 
dimensionless form a spatially constant electrical conductivity was assumed. 

It is true that circulating currents would be impeded by a low conductivity 
near the insulator walls; however, the conductivity does not drop rapidly near 
the cold wall if the plasma has an elevated electron temperature. It was 
shown by Hale and Kerrebrock (ref. 18) that a two-temperature plasma flowing 
over a surface with a normal magnetic field results in a more uniform conduc- 
tivity profile and much larger boundary- layer currents than does an equilib- 
rium flow over the same surface. In some cases the conductivity profiles of 
Hale and Kerrebrock actually exhibit maximum peaks near the wall. 


CONCLUDING REMARKS 


The constant conductivity model of incompressible MHD duct flow predicts 
electrical currents circulating within the working fluid of an open circuit 
J x B device. These currents cause ohmic and viscous heat generation, which 
increase with magnetic field strength and are maximum near the insulator walls 
Heat-transfer rates to the walls of an open circuit J x B plasma accelerator 
were measured; they increased with applied magnetic field strength at the 
insulator walls but remained nearly constant at the electrode walls. This dis 
tribution of heat transfer is qualitatively the same as the calculated heat 
generation distribution. The measured total heat-transfer rates are in good 
quantitative agreement with the calculated total heat generation rates. The 
magnitude and direction of the calculated circulating currents cause signifi- 
cant deceleration of the higher velocity flow near the center of the duct. 
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The problem of the circulating current will exist in any MHD duct in which 
there is an appreciable boundary layer and a large conductivity extending to 
near the walls. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, Feb. 17, 1969 
129-02-08-01-00-21 
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